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Application  of  the  CPEX  Non~Ideal 
Explosive  Model  to  PBXW~115 

DA.  Jones  and  D.L.  Kennedy* 


MRL  Technical  Report 
MRL-TR-91-40 


Abstract 


IVf  present  a  detailed  description  of  the  use  of  the  ICl  non-ideal  detonation  code  CPCX 
to  model  the  performance  of  the  composite  underwater  explosive  P8XW-U5  (20% 

RDX,  43%  AP,  25%  At  and  12%  HTPB).  We  use  both  the  ICl  chemical  eifuilibrium 
code  IDEX  and  Mader's  BKVl  code  to  first  estimate  the  Chapman-Joujfuet  parameters 
for  PBXW-U5.  Two  different  reaction  rate  models  are  assumed  for  the  decomposition 
of  the  explosive  components  and  the  time  constants  fior  each  are  found  by  fitting  to 
experimental  US  data  on  the  variation  of  detonation  velocity  uiJh  charge  diameter.  The 
data  is  best  fitted  by  a  three  term  reaction  rate  model  which  equates  hotspot  content  to 
some  fraction  of  RDX  cpntent,  the  intermedute  scale  reaction  primarily  to  the  AP 
decomposition,  and  the  slowest  process  to  the  Al  reaction.  After  adjusting  some  of  the 
reaction  time  constants  because  of  differences  between  the  US  and  Australian 
compositions  the  model  then  accurately  leproduces  the  variation  of  detonation  velocity  , 
with  charge  diameter  found  for  Australian  PBXW-115. 


*  ICl  Explosives  Croup 
Cate  1,  Ballarat  Road,  Deer  Park.  V'lctorin  J02J 


MATERIALS  RESEARCH  LABORATORY 


I  Accesrtn  tor _ .  , 

I NTIS  CRA4I 


NTIS  CRA&I 

oric  lAB 

U  a-M  Oo  ceJ. 

By..  _ 

Oi  \  ■»  ‘.ic  ! 

\b«i  ,  I  ' 

D*  t  j 

ZJ-JI  i 


•  I 

I 

1 

I 

-1 


2  20  024 


09-04499 

JiEDini- 


PuHisM  hy 

Materials  Research  laboratory 
Cordite  Avenue,  Maribyrnong 
Victoria.  3032,  Australia 

Telephone:  (03)  319  3887 ■  ■ 

Fax:  ,(U3>  318  4536  ■' 

©  Commonwealth  of  Australia  1991 
AR  ^o.  006-833 


AFPROVED  FOR  FVVUC  RELEASE 

'  v 


I 


Authors 


DA.  Jones 

David  lanes  graduated  from  Monash  Unwersity  in  1972 
with  a  BSc  (Hons).  He  obtained  his  PhD  from  Monask 
in  1976.  After  working  at  Strathclyde  University, 
London  University  and  the  University  of  New  South 
Wales  he  joined  MRL  in  1983.  He  has  worked  on  the 
numerical  modelling  of  shaped  charge  warheads  and 
slapper  detonator  devices.  From  February  1987  to  May 
1988  he  IMS  a  Visiting  Scientist  at  the  Laboratory  for 
Computational  Physics  and  Fluid  Dynamics  at  the  Naval 
Res^ch  Laboratories  in  Washington  DC.  While  there 
he  worked  on  advanced  computatiorutl  fluid  dynamics. 


D.L.  Kennedy 


David  Kennedy  graduated  from  Unkersity  of  Melbourne 
with  BSc  (Hons)  in  1976,  and  PhD  in  1979.  He  joined 
Central  Research  Laboratories  of  ICl  Australia  in  1980, 
and  has  been  erigaged  in  long  term  research  into  the  . 
fundamental  delohics  of  the  commercial  explosives 
systems  used  by  mining  and  construction  industries. 
From  1984  to  1987.  he  was  secondedc  to  the  'Cl 
Explosive  Croup  Technical  Centre  in  Scotland.  His 
current  areas  of  interest  include  the  computer  simulation 
of  shock  initiation,  computer  animation,  the  design  of 
initiation  systems  for  blasting,  and  the  hazard  analpis  of 
commercial  explosives. 


Contents 


1.  INTRODUCTION  7 

2.  CPEX  CODE  DESCRIPTION  8 

3.  CALCULATION  OF  EQUILIBRIUM 
STATE  13 

4.  APPLICATION  TO  PBXW-115  16 

5.  DISCUSS'ON  AND  CONCLUSION 

6.  ACKNOWLEDGEMENTS  28 


CHAPMAN-JOUGUET 


2S 


7.  REFERENCES  28 


Application  of  the  CPEX  Non-Ideal 
Explosive  Model  to  PBXW-115 


1.  Introduction 

Several  years  ago  Materials  Research  Laboratory  (MRL)  was  asked  to  advise  on 
the  choice  o<  a  suitable  explosive  tilling  tor  use  in  Australian  naval  mines.  A 
range  ot  explosives  in  current  UK  or  US  service,  or  undergoing  qualification,  for 
underwater  applications  was  surveyed  and  a  recommended  warhead  filling  was 
proposed.  In  addition,  it  was  rerommended  that  R  &  D  work  be  initiated  on 
ammonium  perchlorate  (AP)  based  pwlymer  bonded  explosives  (PBXs)  because 
of  their  reported  lower  vulnerability  and  higher  porformancc.  PBXW-115,  an 
experimental  formulation  uiKler  active'development  in  the  US,  was  identified  as 
bemg  most  suitable  because  of  its  projected  lower  cost  compared  with  other 
PBX  altenutives  and  its  desirable  vulnerability  characteristics.  Another 
attractive  feature  was  that  the  binder  system  was  similar  to  that  used  in 
Australia  for  composite  propellant  production. 

MRL  was  subsequently  tasked  by  Navy  to  consider  the  feasibility  of 
PBXW-IIS  as  a  fill  for  Australian  naval  mines.  Work  on  this  task  to  date  has 
concentrated  mainly  on  formulation  aspects.  PBXW-115  is  a  cast  cured 
explosive  with  a  nominal  composition  of  20%  RDX,  43%  AP,  25%  A1  and 
12%  HTPB  based  polyurethane  binder.  Research  has  concentrated  on  the  blend 
of  RDX  grades,  the  A1  particle  size,  and  both  the  binder  chemistry  and 
mechanical  techniques  for  processing  the  formulation. 

In  parallel  to  these  predominantly  chemical  studies,  we  have  undertaken  an 
investigation  into  the  possibility  of  modelling  the  performance  of  P8XW-115 
and  related  formulations.  This- is  a  new  area  of  endeavour  for  MRL  as 
PBXW-n5  is  a  non-ideal  explosive  ill  and  therefore  not  amenable  to  the  more  ' 
usual  techniques  used  for  the  modelling  of  military  explosives.  Mader.  has 
defined  a  non-ideal  explosive  as  one  having  a  C-)  pressure,  velocity,  or 
expansion  isentrope  significantly  different  from  those  expected  from 
equilibrium,  steady  state  calculations  such  as  those  performed  by  the  computer 
code  BKW  12). .  Johnson  et  al.  (3|  have  modelled  the  performance  properties  of 
a  variety  of  non-ideal  explosives  using  data  obtained  from  the  aquarium  test 


and  have  nude  comparisons  with  results  from  BKW  calculations  assuming' 

.  varying  degrees  of  reaction  at  the  detonation  front.  We  will  consider  this 
approach  in  more  detail  in  a  later  section. 

All  explosives  exhibit  some  form  of  diameter  effect  in  which  the  detonation 
velocity  decreases  with  decreasing  chaige  dianuter  until  a  critica!  dtameter  is 
reached  and  stable  detonation  cannot  be  sustained.  For  ideal  explosives  the* 
change  in  velocity  with  diameter  is  minimal  until  very  close  to  the  critical 
dianuter,  whereas  for  nomideal  explosives  the  diameter  dependence  is  /ery 
pronouitced  and  the  velocity  at  the  critical  diameter  can  be  as  low  as  30%  of  the 
value  at  infinite  diameter.  Commercial  explosives  often  behave  very  non 
ideally  aiul  their  performance  properties  depend  strongly  on  charge  diameter 
and  corrfinement.  These  explosives,  like  military  underwater  explosives, 
contain  separate  fuel  and  oxidizer  species,  often  in  physically  separated  phases, 
and  their  heterogeneous  luture  leads  to  much  larger  reaction  zones  and  more 
curved  detorution  fronts  than  those  in  ideal  explosives  {4]. 

To  model  the  performance  of  conunefcial  explosives  Kirby  and  Leiper  have 
developed  a  snull  divergent  detonation  theory  for  intermolccular  explosives '[5|. 
Their  nxxiel  has  been  developed  into  the  computer  code  known  as  CPEX 
(Commercial  Performaitce  of  Explosives)  and  used  for  a  variety  of  applications 
[6-81. 

In  this  report  we  describe  the  application  of  the  a*EX  code  to  PBXW-115 
using  experimental  data  on  the  detonation  velocity /diameter  effect  obtained  by 
Forbes  et  al.  Ill,  and  also  by  Bocksteiner  et  al.  191  in  a  recent  MRL  field  trial. 
Section  2  describes  the  theory  and  basic  equations  underlying  the  CPEX  code, 
and  Seaion  3  discusses  the  problems  involved  in  calculating  the  equilibrium 
Chapman-Jouguc!  s*’tc  for  PBXW-115,  which  is  treeded  as  input  for  the  CPEX 
code.  In  Section  4  we  describe  the  results  obtained  from  the  application  of 
CPEX  to  PBXW-115,  and  in  the  final  section  we  discuss  these  results  and  our 
approach  in  relation  to  other  methods  which  have  been  used  to  model  the 
behaviour  of  hon-ideal  explosives. 


2.  CPEX  Code  Description 

The  CPEX  code  is  based  on  the  analytical  r)on-idcal  detonation  theory  of  Kirby 
and  Lieper  151,  which  is  a  further  developmenf  of  the  work  of  Own  flOl,  and  is 
ultimately  based  on  the  small  divergent  detonation  theory  of  Wood  and 
Kirkwood  [11].  The  theory  is  applicable  to  both  unconfined  and  confined 
axisyrrunetric  detonations  and  its  application  requires'  the  numerical  solution  of 
a  set  of  coupled,  nbniinear  ordinary  differential  equations  between  the  shock 
front  and  the  end  of  the  reaction  zone. 

A  detonation  wave  is  a  supersonic  shock  w'ave  travelling  through  a  reactive  ■ 
medium.  The  high  temperature  generated  by  the  passage  of  the  shock  through 
the  material  initiates  chemical  reaction  and  energy  release.  A  stable  detonation 
with  a  unique  detonation  velocity  (the  Chapman-jouget  or  C-1  velocity)  will 
then  develop  if  the  chemical  reaction  zone  can  remain  coupled  to  the  wave  and 
continuously  feed  energy  to  the  shock  front.  Any  model  of  a  detonation 
therefore  requires  the  simultaneous  solution  of  the  -".oupled  equations  describing 


both  material  flow  and  chemical  reaction.  In  the  CPEX  model  the  material 
flow  equations  used  are  the  Euler  equations,  a  set  of  three  coupled  partial 
differential  equations  describing  the  conservation  of  mass,  momentum,  and 
energy  for  an  inviscid  fluid.  These  have  the  general  form 

^  ♦  V  '(po)  -  0 


3(pu) 

^  ~  ♦  V'fpot?)  •  -VP 
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where  p  is  the  density,  v  the  fluid  velocity,  and  "  the  pressure.  The  quantity  E 
is  the  total  energy  per  unit  volume  and  is  defined  by 

£  «  pe  ♦  .1  pV’V  W1 


where  E  is  the  specific  internal  energy. 

The  extent  of  chemical  reaction  is  described  by  a  progress  variable  X,  which  is 
the  fraction  of  reacted  explosive  in  the  material.  The  equation  govertting  the 
time  evolution  of  X  has  the  form 


R(p>P,X; 


where  the  exact  functional  form  of  R  will  be  specified  shortly.  Note  that  the 
time  derivative  in  equation  (5)  is  a  total  time  derivative  following  the  fluid  flow 
and  has  the  form  ' 


S’-t  *  (t'-V) 
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The  CPEX  code  tnodels  steady  state  Ui  tonation  in  a  cylindrical  charge  Cor  flte 
case  of  slightly  divergent  flow.  In  this  case  several  simplifications  can  be 
made.  We  first  rewrite  equations  (1)  to  (3)  in  an  axisymme*ric  cylindrical 
coordinate  system 


££  ♦  p(« 

df  ^  * 


♦  w.)  ♦  p  _«.  •  0 
:  r 


(7) 


(8) 


*  P. 


-  0 


(9) 


d«  P  do 

Tt  '7  *di' 


(10) 


!iere  u  and  w  are  the  axial  (t)  and  radiai  (r)  components  of  the  velocity  vector  v 
and  SLjscript  *  or  r  denotes  partia  if/erentiation  with  respect  to  i  or  r. 

Using  the  steady  state  condition  and  the  assumption  of  small  divergent  flow, 
and  considering  motion  only  in  the  axial  direction  (at  r  =  0,  i.e.  along  the 
centreline),  equations  (7)  to  (1C)  and  equation  (5)  can  be  written  in  the  form 
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Kirby  and  Lieper  [5]  further  simplify  the  above  set  of  equations  by  using  the 
equation  of  state  in  the  form  e  -  e(p,  P,  \)  to  derive  the  following  expression  for 

K 


where  c  is  the  local  sound  speed  and  o  is  die  thcrmicity,  defined  b>' 
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the  quantity  o  is  a  measure  of  the  heat  release.  The  systtm  is  exothermic  if « 
is  positive  and  endothermic  if  <t  is  negative.  Use  of  equation  (16)  then  results 
in  the  following  definii.g  set  of  equations 
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The  term  in  equations  (18)  and  (19)  represents  the  radial  divergence  of  the 
flow.  If  the  flow  was  completely  one-dimensional  this  term  would  be  zero  and 
the  equations  would  represent  the  standard  Zeldovich,  von  Neumann,  Doering  ■ 
(ZND)  model  of  detoriation  (12).  In  this  model  the  initial  disturbance  shocks 
the  material  to  the  von  Neumann  point,  which  results  in  an  increased  value  of 
the  density,  and  considerably  higher  values  of  both  pressure  and  temperature. 
This  irutiates  chenucal  reaction,  which  then  proceeds  at  a  finite  rate  until  all  iif 
the  explosive  has  been  converted  to  detonation  products.  In  the  qne- 
dimeruional  case  (i.e.  w,  =  0),  equ4ticns  (18)  to  (21)  can  be  used  to  show  that 
the  end  of  the  reaction  zone  always  coincides  with  the  Chapman-Jouget  (QJ) 
point,  which  is  defined  as  the  point  at  which  the 'sum  of  the  fluid  velocity  plus 
the  local  sound  speed  is  equal  to  the  detonation  velocity.  This  is  an  important 
point  in  the  flow  because  if  reaction  occurs  beyond  this  point  then  the  energy 
which  is  released  is  unabiO  to  coiitribute  to  the  support  of  the  detonation  wave. 

For  divergent  flow  the  simple  ZND  picture  no  longer  applies,  and  some 
degree  of  reaction  will  always  occur  beyond  the  C]  point.  Equations  (18)  to 
(21)  (for  finite  w,)  can  be  used  to  show  that  the  reactjon  must  satisfy  the 
following  criteria  at  the  CJ  point 

•  iw./u  (22) 

T 

which  implies  that  the  rate  of  chemical  energy  release  is  balanced  by  the  energy 
diminution  to  ‘he  lateral  flow  expansion.  Thus  for  radially  divergent  flow 
there  are  two  sources  of  energy  loss:  (a)  the  energy  loss  to  the  radial  motion, 

(b>  the  fraction  of  chemical  energy  released  beyond  the  C/  point.  This  means 
that  the  detona.ioi'  velocity  for  a  cylindrical  charge  with  a  finite  diameter  will 
always  be  less  than  th."*  planar  or  CJ  value,  and  will  only  approach  the  CJ  value 
as  the  diarr  eter  becomes  infinite. 

Before  equations  (18)  to  (21)  can  be  solved  an  expression  for  the  radial 
divergerKe  w,  must  be  specified.  It  is  easily  shown  (121  that  le,  can  be  related 
.  to  the  radius  of  curvature  of  the  detorratidn  front,  R^,  by  the  expression 

Wj  •  (D  *  u(o))/Rj  (23) 

where  D  is  the  detonation  velocity  and  iiio)  is  the  flow  velocity  at  the  shock 
front.  The  problem  thus  reduces  to  the  specification  of  Rj.  The  developers  of 
the  CPEX  code  surveyed  maiiy  Explosives  and  found  that  Rj  was  irtversely 
proportional  to  charge  diameter  d,  provided  that  both  R^  and  4  wne  srai^  by 


projx)rticnality  constant  depends  on  whether  the  explosive  is  a  Group  1  or 
Croup  2  explosive  (as  defined  by' Price  (13)),  but  the  correlation  is  sufficiently 
Strong  enc’igK  to  provide  the  r  3cessary  relationship  between  Rj  and  d.  If  more 
e>.act  experimental  data  are  available  for  a  particular  explosive  then  CPEX  has 
an  option  to  iccept  this  data.  We  used  both  methods  in  our  analysis  of 
PBXW-115. 

The  reaction  rate  law  used  by  C!PHX  has  the  following  functional  foin  to 
allow  for  the  multi-phase,  multi- component  nature  of  composite  explosives 


dX 
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The  subscripts  h,  I,  and  s  refer  to  the  hot  spot,  matrix  (in  commen.'^  explosives 
this  is  often  a  liquid  phase)  and  included  solid  phases  within  the  explosive- 
The  pressure  expone''ts  (i  =  h,  I,  s)  are  normally  set  to  unity  and  -the  flj  are 
gaussian  form  fijnetions  which  switch  the  reactions  on  and  off  as  the  various 
phases  are  ignited  and  consurned.  The  depend  on  X  and  the  initial 
formulation  of  the  explosive.  The  tj  are  characteristic  reaction  times  for  each  of 
the  phases,  and  P^^  is  «.  hot  spot  critical  pressure.  In  typical  applications  of 
CPEX  the  dj  and  bj  are  fixed  and  the  T|  and  P),  are  ustxl  as  adjustable  fitting 
parameters. 

The  equation  of  state  for  a  single  phase  has  a  simple  density  dependent 
polytropic  equation  of  the  form 


«  -  £(g-l)-> 
P 


(25) 


g  ^  So  *  StP  *  S2P^ 

The  reasoning  behind  the  choice  of  equatiens  (25)  and  (26)  to  desbibe  the 
equation  of  state  is  discussed  in  detail  by  Kirby  and  Leiper  f5’].  The  values  of 
the  constants  tor  the  explosive  products  are  found  by  fitting  to  isentrope  data 
from  an  ideal  thermodynamic  code,  while  the  constants  for  .he  unreaefed  phase 
are  obtained  by  considering  the  sho,ck  state  of  the  explosive  a;.d  using  known 
Hugoniot  data.  Tlie  unreKted  solid  state  is  usually  a  mixture  of  ingredients, 
and  in  this  case  the  Hugoniots  of  the  components  arc  combined  using  the 
riiethod  of  >fanasenkov  et  al.  (141. 

Equations  (18)  to  (21)  are  a  coupled  set  of  nonlinear  ordinary  differential 
equations  which  must  be  solved  subject  to  a  mixed  set  of  boundary  conditions; 
at  the  CJ  plane  both  the  CJ  condition  and  equation  (22)  must  be  satisfied,  and 
immediately  behind  ,thc  shock  front  the  variables  mu^  agree  with  the  uhreacted 
shock  state.  The  detonation  velocity  now  becomes  an  eigenvalue  of  the 
equation  set  and  its  value  is  dependent  on  the  charge  diameter,  dec  ree  of 
.  cordinement,  and  on  fte  chemical  reaction  rate  of  the  explosive. 

In  typical  applications  of  CPEX  the  equations  o"  state  of  both  reactants  an-J 
products  are  known  artd  the  ideal  C)  parameters  are  obtained  from  a  chemical 


equilibrium  code  such  as  BKW  or  the  IQ  IDEX  code.  The  reaction  rate 
parameters  T,  and  P/,  are  then  estimated  by  using  CPEX  to  obtain  a  fit  to 
experimental  data  on  the  variation  of  detonation  velocity  with  ch^e  diameter, 
and  then  CPEX  can  be  used  to  predict  the  explosive  performance  for  a  variety 
confinements  [7],  as  well  as  the  effect  of  chaige  diameter  and  porosity  on 
ini'^ation  behaviour  [6],  and  also  to  make  comparisons  witli  other  theoretical 
approaches  for  predicting  the  .perfor.nance  of  condensed  explosives  [8} 


3.  Calculation  of  Equilibrium 
Chapman-Jouguet  State 


Before  running  CPEX  on  a  particular  explosive  formulabon  an  estimate  must  be 
made  of  the  CJ  velocity  and  ideal  CJ  gamma  for  the  given  explosive.  This  can 
be  done  using  any  chemical  equilibrium  detonation  code  and  we  have  used 
both  Mader's  BKW  code  and  the  ICI  IDEX  code  to  estimate  these  parameters  for 
PBXW-;i:. 

The  BKW  code  employs  the  Becker-Kisliakowsky-Wilson  (BKW)  equation  of 
state  to  calculate  the  detonation  properties  of  explosives  [15].  Steady-state 
time-independent  behaviou’  is  assumed  and  the  detonation  products  are 
considered  to  be  in  instantaneous  chemical  equilibrium.  For  an  explosive 
consisting  of  a  mixture  of  chemical  com|X>unds  the  necessary  input  data  are  the 
fraction  by  weight,  elemental  composition,  and  heat  of  formation  of  each 
compound.  The  code  requires  the  entropy  of  the  product  species  to  be 
cxpicssed  as  a  polynomial  function  of  temperature  and  th.n  the  detonation 
parameters  are  calculated  by  iteratively  minimizing  the  free  energy  once  an 
initial  guess  of  the  product  composition  has  been  made. 

The  MRL  version  of  BKW  is  run  through  a. user  interface  package  known  as 
'JjcRBKW  [16].  'niis  is  an  interactive  FORTRAN  program  which  prepares  an 
input  data  file  for  BKW  by  asking  a  series  of  questions  about  the  explosive  and 
its  products.  USERBKW  also  maintains  three  database  ftl«,  BKWCW, , 
BKWGAS  and  BKWkDL.  BKWCMP  contains  information  on  standard 
explosive  components  while  BKWGAS  and  BWKSOL  contain  information  on 
gaseous  and  solid  detonation  products.  It  was  rtecessary  to  uprfate  BKWCMP 
by  the  addition  of  the  elemental  composition,  heat  of  formulation,  dettsity,  and 
formula  weight  for  AP,  while  BKWGAS  required  the  heat  of  formation,  entropy 
constants^  and  the  covolumes  for  both  HQ  and  Q2.  These  were  obtairred  from 
standard  compilations  {17].  For  the  purposes  of  tire  BKW  calculation  the  MTPB 
binder  in  reXW-115  wa«  replaced  by  an  existing  plasticizer  in  die  BKWCMP 
database  This  should  have  negligible  effect  on  me  calculated  Q 

parameters. 

USERBKW  provides  three  methods  for  estimating  the  initial  amount  of  each 
of  the  assumed  detonatior  products;  explicity,  where  the  user  gives  his  own 
estimates,  fiierarchically,  where  the  user  assigns  priorities  with  first  priority 
being  pven  fo  the  proc.  ct  most  Ikely  to  be  formed,  and  by  default  We  used 
both  the  erlplicit  and  hi-Tarchical  methods  and  these  gave  quite  different  initial 
estimates  of  the  equilibrium  prodiicts.  However,  when  both  these  estimates 
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niTC  jjiv-'n  to  BKW  and  tho  (rt’O  energy  minsmi/ation  procedure  was  followed 
ihe  final  ecjuilihr'.um  state  prixiuced  by  each  estimate  was  identical,  as  Table  1 
•.hoiv-  :  The  detonation  parameters  at  a  density  of  1  80  Mg/m’  calculated  by 
HisTV  are  also  listed  in  lable  1. 


Zjhlf  1:  D(t.y.utu  r.  rrejacts  'er  PBXVV1I5  CjuuUuJ  by  BKW  Code 
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,  The  C'l  pressure  of  TH  *i  (.Pa  and  deton.ilion  velocity  of  8  Ot  mm/ps  are  both 
biftber  .lun  esjxsled  ortxs  Kis  nx'jsund  llw  (  !  pnssore  of  PbXW  Il^ 
using  the  aijuaniim  lisi  nx  tfxd  and  fourvf  a  valm*  clow  to  1?  2  C.I’a.  While 
the  anaivsis  of  tfk-  ainiarmm  tist  data  is  bas«d  on  espnssions  which  are  only 
valid  for  idi’al  esplmms,  s<i  ih.it  some  ermr  is  I'nindund  in  calculating  (T 
pnssuns  for  m  n  ideal  e«plosi\>s.  n  is  thought  that  this  would  inhiHlucc  errofsi 
of  no  more  than  aNmt  ,  vi  die  value  of  2H/H  (.Pa  must  be  viewed  with  v'me 
suspmon 

I  oflxs  has  also  rtv.isun  d  the  detx'ixfense  of  lt«'  deforvifion  velocity  on 
I  Kirge  diameter  1 1 1,  an»l  Ihi’. •■•if.ipolalid  infinite  ifi  imeter  detonation  veiiwitv 
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from  this  data  is  6.195  nun/ps;  the  value  obtained  from  the  .vfRL  detonation 
velocity /diameter  measurements  is  5.650  mm/ps.  These  values  may 
und^restintate  the  CJ  velocity  of  PBXW-115  because  in  both  cases  the  maximum 
charge  radius  used  vvas  only  approximately  twice  the  critical  diameter,  and  this 
may  have  prevented  some  late  bme  reactions  from  contributing.  When 
considered  in  conjunction  with  the  discrepancy  between  the  measuntd  and 
predicted  CJ  pressure  however,  the  difterence  between  the  extrapolated  and 
predicted  C)  velocity  also  casts  some  doubt  on  the  validity  of  the  GJ  slate 
calculated  by  BKW, 

We  have  also  esbmated  the  CJ  state  of  PBXW-tl5  using  the  IGI IDEX 
chemical  equilibrium  code  and  the  results  are  shown  in  Table  2.  luEX  uses  an 
intermolecular  potential  equation  of  sute  for  fluids  and  a  Mumaghan  equation 
of  state  for  solids. 


Table  2:  Detonation  Products  for  PBXW-17S  Calculated  by  IDEX  Code 


Species 

Mole  Numbers 
(100  g  esplosise) 

CH, 

0.1556 

CO 

0.5436 

CO3 

000673 

H: 

027T5 

.VH, 

0  4661 

H;0 

02681 

HO 

0  2591 

N'l 

02197 

NO 

0A5.W  -  KT* 

0- 

0  2524  .  UT' 

ci  . 

0C5340 

C  (Graphite) 

1.0  «  10’' 

C  (Diamond) 

06338 

A1;0, 

0  4632 

Al 

I  0  •  10  -’ 

C)  Parameters: 

Delonajioo  Pn-ssure  «  22  5!SCPa 
Detonation  Velocity  =  teeiOiil  m/s 
Di'lonatlon  Tcmperalurc  •  52V7  K 
Density  ■  t.7W)Og/cm’ 

Gamma  •  ■  2  5268 


The  calculated  CJ  pressuir  of  22 .5  GPa  is  still  consKierabiy  higher  than  the 
exjii'nmental  value,  but  also  significantly  kiwer  than  the  value  calculated  by 
HKW,  while  the  Cl  svV'Kity  of  ft  ft6  mm/ps  is  tiow  much  cIomt  to  the  value 
obtained  frt>m  the  exirapi>latK>n'of  the  expenmenlal  detonation  ■ 
velooly/diametiT  data. 


A  comparison  of  the  equilibrium  compositions  predicted  by  the  two  codes 
shows  considerable  differences,  with  the  only  clear  point  of  agreement  being 
complete  conversion  of  all  A1  into  0.463  moles  of  AUOj  per  100  g  of  explosive. 
Whilp  both  codes  show  appreciable  amounts  of  carton  present  in  the 
equilibrium  composition,  BKVV  has  predicted  that  this  will  be  in  the  graphite 
form,  while  IDEX  has  predicted  that  most  carbon  will  be  present  in  the 
diamond  form. 

We  have  used  both  estimates  of  the  CJ  state  as  input  to  the  CPEX  code 
calculations  described  in  the  next  section. 


4.  Application  to  PBXW~115 

As  dcscnbcd  in  Section  2,  the  reaction  rate  model  used  in  CPEX  involves  three 
dist-nct  stages,  each  with  a  characteristic  reaction  time.  ■  In  the  initial  ignition 
phase  the  shoe*  front  compresses  the  explosive  and  produces  hot  spots.  The 
ho*,  ^pots  then  initiate  a  fast  reaction  in  the  matrix  (or  liquid  phase),  and  this  is 
followed  by  a  slower  reaction  of  the  solid  phase  as  the  matrix  is  depicted. 

CPEX  was  developed  for  a’mmerrial  explosives  and  these  typically  contain  a 
large  friction  of  voidage.  of  the  order  of  15"  to  30'«  in  many  cases.  These 
voids  are  important  as  iney  are  considered  to  bt>  the  sites  at  which  the  hot  spots 
are  formed  as  the  inmaimg  shock  compresses  the  explosive.  In  PB.\W-115, 
however,  great  care  is  exerased  in  processing  to  ensure  that  the  voidage  is 
minimi7c>d.  Compositions  recently  produced  at  MRl.  for  example  contain 
yoidage  levels  of  less  than  2" . 

In  order  to  apply  CPEX  to  PBXVV-1I5  it  is  first  mvessary  to  identify  each  of 
the  three  .stages  in  the  reaction  rale  model  with  a  corresponding  process  in  the 
explosive  decoirposition.  Two  possibilities  are  immediately  apparent,  one  is 
to  asenbe  the  hot  spot  formation  to  the  v<>idage.  as  in  commt'rcial  explosives, 
and  then  to  attribute  the  fast  reaction  to  ifx-  RDX  decomposition  and  the  slower 
reaction  to  the  A  P  and  A'  decompi>sition.  The  other  possibility  is  to  assume 
negligible  voidage  and  equate  the  hot  spot  content  to  some  fraction  of  the  RDX 
cimlenl,  the  intermediate  scale  reaction  primarily  to  the  AP  decomposition  and 
fhe  slowest  process  to  the  Al  reaction.  We,  pcTfonrwxJ  calculations  for  both 
schemes  and  the  results  are  describi'd  in  this  section 

1  he  first  results  we  describe  are  those  in  which  we  considcTed  the  RDX 
component  to  correspond  toi  the  fast  reacting  liquid  phase,  while  the  AP,  ,A1  and 
binder  together  constilulcsd  the  slower  burning  solid  phase.  We  used  a  density 
of  1  80  Mg/m'  for  RDX  and  1  83  Mg/m*  for  the  combined  AP,.Arand  binder; 
the  PBXW-115  had  a  density  of  I  78  Mg/m*  (this  is  the  density  of  the 
PBXW-115  used  by  Forties  rf  a/  (11).  This  ltd  to  approximately  6%  voidage, 
which  is  certainly  higher  llvin  would  occur  in  practice.  These  calculations  are  . 
not  intended  to  pnduce  exact  numbers,-  they  are  ustd  pnmanly  to  illustrate, 
the  capabilities  of  the  CPEX  code  and  explore  its  applicability  to  PBXW-115  and 
similar.  undtTwater  explosives,  faweral  of  the  inpqt  params’tcrs  are  specified 
(inly  approximately,  but  with  soIfKwnt  precisiois  so  that  confidcmce  can  be 
p'actd  in  the  gtmeral  trends  shown  by  the  results. 


The  input  information  required  to  run  CPEX  is  shown  in  Table  3.  The  ideal 
detonation  velocity  and  CJ  gamma  were  obtained  from  the  BKW  results 
described  in  the  previous  section,  and  the  heat  of  reaction  was  estimated  from 
the  standard  relationship  between  detonation  velocity,  gamma,  and  heat  of 
reaction.  The  Hugoniot  parameters  for  RDX  were  taken  from  the  Lawrence 
Livermore  compilation  (18|,. while  the  Hugoniiil  parameters  for  the  solid  phase 
were  combined  using  the  method  of  Afanasenkov  et  at  |I4|  from  the  parameters 
for  AP  1191,  A1  1301  and  binder  1211.  The  hot  spot  reaction  time  and  criticai 
pressure  and  liquid  and  solid  reaction  times  were  estimated  by  fitting  to  the 
experimental  data  of  Forbes  ei  ul  dc'senbing  the  vanation  of  detonation  velocity 
with  charge  diameter  for  unconnned  charges. 

Table  i:  Infut  Data  t-r  CPEX  Ci'Je:  hot  ipot  content  eaiaiteJ  to  initial  ixijage.  Cl 
parameters  from  fiklV  ci  Je.  anJ  Jctonalwn  :vl,xtty  mcasitrcmentf  fr^tnt  Forbes  et  jl. 
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Figure  I  shows  the  detonatKin  veliKity  •  inverse  diattmer  relationship  for 
unconfined  PBXW-1I5  predicted  by  CPEX  using  the  pararrH-ltT  values  in 
Table  3.  The  lower  branch  of  the-  curve  represents  an  unstable  set  of  scilutKins 
to  the  divergent  flow  theory.  It  coiresponds  to  sets  of  omdihon*  which  can 
just  maintain  a  steady  detonation  in  the  absence  of  perturbations,  but  which 
wiiuld  eithcT  dtvay  to  failure  or  U-ad  to  detonation  at  tnc  high  vilue  in  the 
presence  of  any  random  fluctualKins.  The  uppvT  branch  of  the  ''irVc  is  a 
reasonable  fit  to  the  expi-nnx>ntal  data.  The  pndicted  critical  diai.w-tcr  is 
34  5  mm,  which  is  vcTy  cliwe  to  the  expi’nnwntal  value  of  37.1  »  1.6  mni  !t|. 
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fiifurr  I.  {V!,  itUoiv  inirTV  char j(e  Jutmriir  fvr  uncon  fine  J 

I'HSW  115  THc  (rurs^-.'.-s  a'f  Jala  from  for^YS  et  al  III.  the  sold  line  ii 

th(-  CPfX  •:!  uMr<  ;v  Jjlj  m  Tutle  3. 

(orbrt  rt  jI  hjM-  al«o  mcasurvd  the  dotoiutKm  velotity  •  durrftrr  i-ffcct  for 
ih.ir>;i->  fonlinixl  in  2  5  mm  Jind  5  0  thm  bras*  tubes.  They  found  that  the 
intK  jl  dtameter  fiir  '  mm  thicK  bras*  confined  charj^es  wa*  about  one  half  that 
ot  the  unconfin»Nl  chtica'  '.ameter.  ,  Since  the  difference  between  the  enheal 
dumeters  of  2  5  and  5  0  mm  bras*  confined  chaises  was  *o  small  they  expected 
thit  only  slightly  smaller  critical  diameters  would  be  obtained  if  brass 
confinc'ment  greater  than  5  0  mm  was  used.  Figure  2  shows  the  CTEX 
predictKm  for  the  deto^tion  veloaty  -  diameter  effect  for  charge*  confined  iii 
'  bras*  tubes.  The  CTEX  calcutatiofi  assumes  an  effectively  infinite  degrei*  of 
confinement,  and  approximate  shock  transit  time  calculations  show  that  5  mm' 
thick  brass  confin^Tnent  provide*  sufficient  impedance  to  prevent  the 
rarefactions  fnim  effecting  the  reaction  zone.  The  agreement  is  not  quite  as 
giHHl  in  this  caw:  CTEX  predict*  a  critical  diameter  of  22  mm,  while  the 
experimental  value  is  approximately  17  mm. 


As  discussed  in  Section  2,  when  considering  two-dimensionjl  divergent 
detonations  some  concepts  from  the  one-dimensional  hydrodynamic  theom  of 
detonation  require  nvidification.  One  of  these  is  the  concept  of  the  CJ  point. 

In  one-dimensional  theory  the  CJ  point  is  defined  as  the  point  in  the  flow  at 
which  the  sum  of  the  particle  velocity  and  local  sound  velocity  is  equal  to  the 
detonation  velocity.  At  this  point  the  rate  of  reaction  is  zero,  and  the  reaction 
zone  can  be  defined  as  the  distance  between  the  shock  front  and  the  Cl  point 
In  two-dimensional  flow  the  CJ  point  must  satisfy  two  conditions;  the  first  is 
the  one-dimensional  condition  defined  above,  the  second  is  due  to  the  , 
divergence  of  the  flow  and  requires  that  the  reaction  satisfy  a  criterion  such  that 
the  rate  of  cher  .ical  energy  release  is  balanaxl  by  the  energy  diminutfon  due  to 
lateral  flow  expansion  (i.e.  equation  (22)).  Thus  the  reaction  rate  will  alwavs 
have  a  finite  value  at  the  CJ  point  in  divergent  flow  and  the  identification  of  the 
reaction  zone  will  be  less  precise.  This  eftect  is  illustrated  in  Figure  3,  which 
shows  the  extent  of  reaction  at  the  CJ  point  as  a  funclum  of  the  inverse 
diameter.  At  a  diameter  of  41  mm  ft>r  example  this  plot  shows  that  only  23  ' 
reaction  has  txcurred  by  the  C|  point.  As  the  diameter  increases  the  flow 
bevomes  rrusri  oiv-dimensional  and  we  find  the  extent  of  riMction  tending  to  Hi 
at  the  CJ  pt'int  as  the  diameter  bexromes  infinitely  large. 
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rfjprre  2;  Detonation  tirliV'fy  ixtsus  rntrrsc  charge  Jumctiv  for  fh’.VlV  II.s  con'ireJ 
in  bras*.  The  squares  are  ezjvhmen'ai  data  from  Torbes  el  at  HI.  the  se/iJ  Jmc  is  the 
CI’EX  fit  using  the  data  in  Table  3. 


3. CIS  C.O.I’O  0.0:5  0.030 


:h3r;e  Sia.neter. 

Figure  J;  Cxtfr.t  c'  roicHcn  at  the  Cl  point  prrsu..  inverse  charge  diameter  for 
w.i.  r.nncJ  PBXW  US  calculated  by  CPEX  using  the  input  data  in  Table  3. 

I  ho  Cl  Uistatico  IS  dofiiH'd  as  the  di'.tanec  between  the  shock  front  and  the  CJ 
p'int  and  rij;ure  4  >ht>w  s  a  plot  of  Cf  distance  against  inverse  diameter.  At 
infinite  diameter  this  is  etqual  to  the  ide-l  one-dimensibnal  reaction  zone  length 
and  Figure  4  indicali-s  that  this  is  appro'iipately  70.  mm.  As  the  diameter 
decreases  the  C|  distance  dosTeases  as  well  because  more  of  the  reaction  energy. 
IS  contributing  to  dnving  the  '  itf-ral  divergctKC.  Figure  4  indicates  CJ 
distances  pf  approximately  !  nm  to  13  mm  in  the  diameter  range  40  mm  to 
70  mm.  Forbes  et  al  hase  usesj  .measuieirients  of  the  radius  of  culrvaturc  01  the 
detonation  front  and  used  the  original  Wood-Kirkwood  theory  to  estimate 
reaction  zone  lengths  of  2  to  4  mm  for  this  diameter  range.  The  latter 
estimate's' of  CJ  distances  differ  from  those  of  CFEX  due  to  the  effects  of  several 
assumptums  made  in  the  original  Wood-Kjrkwood  theory.  These  included'a 
simpllfu'd  geometry  of  the  reaction  zone,  an  artificial  pressure  profile,  density 
ratios  behind  the  shock  and  al  the  sonic  plane  appropriate  to  a  particular  liquid 
explosive,  and  the  assumption  that  reaction  is  complete  at  the  sonic  plane.  The 
II. finite  diameter  reaction  zone  length  61 70  mm  is  large  mainly  because  the 
roaclinn  is  much,  slower  than  those  fur  the  common  ideal  military  explosives. 
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Figure  4:  Distance  between  shock  front  and  CJ  point  tersus  innrse  charge  dumctcr 
for  unconfined  PBXW-U5  calculated  by  CPEX  using  the  input  data  in  Table  .V 

■  Because  of  the  very  low  level  of  voidage  in  PBXW-115  a  more  likely  approach 
is  to  assume  that  hot  spots  fonn  at  RDX  sites  in  the  matrix,  and  thus  to  equate 
■  the  hot  sfiots  term  in  the  reaction  rate  law  with  some  fraction  of  the  RDX 
content  of  the  explosive.  We  have  found  improved  agreement  with,  the  NSWC 
data  by  considering  a  model  in  which  half  of  the  RDX  is  used  to  form  hot  spots 
and  the  remaining  RDX  and  all  of  the  AP  ^nd  binder  react  on  the  mtermediatc 
time  scale,  with  the  Al  alone  reacting  on  the  slowest  tinte  scale.  VAk.>.have  also 
used  the  CJ  parameters  obtained  from  the  IDEX  code  for  this  particular 
calculation,  and  a  complete  list  of  the  input  data  for  the  CPEX  run  is  shown  in 
Table  4. 


Table  4:  Input  data  fhr  CPEX  Code:  hot  spot  fraction  equated  to  50%  initial  RDX 
content.  Cl  parameters  from  ICI  IDEX  Code.  Detonation  velocity  measurements  from 
Forbes  et  al.  III. 


Explosive  porous  density 

g/cm* 

1  ';900 

Ideal  detonation  velocity 

km/s 

to350 

Heat  of  rcactiuh 

MJ/kg 

8.2(XX) 

Ideal  Effective  Energy 

MJ/kg 

4.6740 

Gamma  Ideal  ^  expanded  gas 

IJ4I0 

Gamma  Cj  •  ideal  CJ 

16490 

Density  liquid  ingredients 

g/cm^ 

1.6870 

Hugoniot  intercept  -  liquid 

km/s 

2.4610 

Hugeniot  slope  -  liquid 

1.6000 

Density  solid  ingredients 

g/cm* 

2.70a’ 

Hugoniot  intercept  •  sohd 

km/s 

5.3500 

Hugoniot  slope  -  solid 

1.3400 

Mass  fraction  of  liquid 

0.7500 

Hot  ‘‘pot  reaction  time 

MS 

.  34.2C(» 

Hot  spot  critical  pa*ssure 

UPa 

,  3J000 

Liquid  reaction  time 

MS 

380,1XX)0 

S>!id  reaction  time 

MS 

lOiUWO 

Liquid  hot  spot  factor 

o.uw 

S>lid  hot  spot  factor 

iroixio 

Kinetic  model  flag  (default  s  (’) 

1.0(/00 

Hot  spot  pressure  exponent 

2.0030. 

L.quid  pa-ssure  exponent 

1.0000 

S4)lid  p.essure  exponent , 

1.0000 

.  Fif;urc  5  -.hows  the  prodictiij  dcKmation  velocity  as  a  function  of  inverse 
diameter  for  the  input  data  in  Table  4.  This  shows  a  much  better  fit  to  the 
experimental  data  (cf.  Fig.  1)  and  predicts  a  critical  diameter  of  38.5  mm,  which 
IS  within  the  experimental  value  of  37.1  ♦  1.6  mm.  Figure  6  shows  CJ  distance 
versus  inverse  diameter. 

Bot  ksteiner  et  al.  [91  have  recently  performed  a  series  of  detonation  velocity 
measumments  on  cylinders  of  Australian  PBXW-115  with  diameters  in  the 
range  40  mm  to  200  mm.  The  results  (along  with  a  CPEX  fit  to  the  data)  arc 
"shown  in  Figure  7.  The  critical  value  for  the  Australian  made  explosive  was 
found  to  be  approximately  80  mm,  which  is  about  twice  the  value  for  the'US 
composition.  Hence  both  the  US  and  Australian  measurements  covered  a 
range  of  approximately  one  to  two  critical  diameters  for  their  respective 
formulations,  and  within  this  range  the  detonation  velocities  are  very  similar. 
Figun.'  8  shows  both  sets  of  data  as  a  function  of  diameter  d,  when;  d  has  been 
divided  by  the  appropriate  critical  diameter  d„  for  each  explosive.  For  a  given 
value  of  scaled  diameter  the  two  sets  of  results  differ  by  no  more  than  a  few. 
percent. 
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Figure  5:  Detonation  ivlocity  oersus  inverse  charge  JiameUr  for  un(Ofi‘ir.o.l 
PBXiV  US.  The  triangles  are  the  experimental  data  pom  Forivs  et  al.  III.  the  solid 
line  is  the  CPEX  ft  using  the  data  in  Table  4. 
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Inverse  charge  diameter,  i/.mm 

Figure  6;  Distance  between  shock  poni  and  Cl  point  versus  inverse  charge  diametv- 
for  unconpned  PBXW-U5  calculutei  by  CPEX  using  the  input  data  in  Table  4. 
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figure  7:  DitiV’.atk'n  tvh’city  ver^uf  inverse  charge  diameter  for  unconfineJ 
rf'  V'.V  115  The  triar.eh'S  are  the  ex^hrimental  data  from  Bxksteiner  et  al.  191,  the 
>olU  line  is  the  CPTX  ‘it  using  the  data  in  Table  4,  and  thr  changes  described  in  the 
text. 


While  the  Australian  and  American  formulations  arc  nominally  the  same, 
there  an?  differences  in  the  nature  of  the  RDX  used  in  the'  two  compositions 
because  of  differences  in  the  method  of  manufacture.  The  American  RDX 
contains  approximately  8%  HMX,  while  the  Australian  RDX  is  completely  HMX 
frc’e.  The  HMX  content  would  nrtake  the  US  PBXW-llS  slightly  more  shock 
sensitive,  and  hence  may  explain  the  smaller  critical  diameter.  The  different 
sensitivities  of  the  RDXs  used  in  the  two  compositions  will  also,  change  tfte  time 
constants  in  the  reaction  rate  law  for  those  processes  connected  with  the  RDX 
decomposition.  'To  obtain  the  CPEX  fit  to  the  data  shown  in  Figure  7  we 
increased  the  hot  spot  reaction  time  from'34  ps  to  110  ps  and  the  intermediaie 
reaction  time  from  380  ps  to  800  ps.  The  valtKis  of  all  other  constants  were  the 
same  as  used  in  obtaining  the  CPEX  fit  to  the  NSWC  data  shown  in  Figure  5. 
The  excellent  fit  to  the  data  shown  in  Figure  7  is  ertcouraging  evidence  that  we 
arc  correctly  modelling  the  basic  rciction  mechanisms  omrring  in  PBXW-llS. 
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Figure  8:  Detonation  velocity  as  a  function  of  scaled  inverse  diameter  for  the  data  of 
Forbes  et  al.  Ill  (triangles)  and  Bocksteiner.  et  al.  191  (squares). 


Wo  intend  to  further  validate  this  model  by  comparii  g  its  predictions  against 
a  variety  of  experimental  data.  The  effect  of  tx)th  brass  and  steel  confinement 
on  failure’  diameter  will  be  calculated  and  compared  with  the  results  of 
experiments  currently  in  progress  at  MRL  (Bockstciner  (221),  and  we  also  intend 
to  use  a  version  of  the  DYNA2D  code  containing  the  CPEX  reaction  rate  scheme 
and  equation  of  state  to  model  aquarium  test  data  obtainea  using  the  US 
formulation.  A  reactive  Flux-Corrected  Transport  code  will  then  be  used  to 
simulate  the  comer  himing  experiments  of  Forbes  et  al.  (I]. 


5.  Discussion  and  Conclusion 

I 

Tne  application  of  CPEX  to  PBXW-llS  enables  us  to  calculate  many  explosive 
properties  as  a  function  of  the  charge  radius.  It  also  allows  us  to  obtain 
information  on  the  reaction  time  scales  of  the  various'  cornponents  and  the 
degree  of  reaction  at, the  CJ  point.  CPEX  provides  a  complete  model  for  a 
slightly  divergent  non-ideal  detonation  for  explosives  of  commctcial  interest. ' 
Other  methods  have  been  published  which  estimate  the  performance  of 
commercial  or  non-ideal  explosives  and  it  is  pertinent  to  mention  some  of  these 
here. 

Johnson  et  al.  (31  have  used  data  from  the  aquarium  test  and  combined  this 
with  BKW  and  two-dimensiorul  code  calculations  to  infer  infbrmalioi:  about  the 
degme  of  chemical  reaction  at  the  detonation  front  for  several  commercial 
explosives.  A  BKW  calculation  was  made  using  the  iiritiai  de' c  ity  and 
chemical  composition  of  the  explosive  assuming  thermodynamic  xctuilibrium 


and  complete  chemical  reaction  at  the  detonation  front.  The  BKW  calctJations 
were  then  compared  with  the  measured  detonation  velocity.  Generally,  the 
measured  detonation  velocity  was  welt  below  the  theoretical  prediction, 
indicating  that  reaction  was  incomplete  at  (or  vety  near  to)  .'-a  detonation  front. 
Additional  BKW  calculations  were  then  performed  with  varying  amounts  of 
the  constituents  withheld  from  chemical  reaction  at  the  detonation  front  until 
agreement  was  achieved  between  measured  and  calculated  detonation  speed. 

If  the  explosive  contains  several  non-energetic  components  then  some  chemical 
knowledge  is  required  at  this  stage  to  decide  whi.h  components  should  be 
reacted.  An  equation  of  state  of  the  detonation  pnxlucts  can  then  be 
constructed  and  used  in  a  two*dimensional  hydrocode  calculation  to  simulate 
the  aquarium  test  data  and  predict  the  position  of  measured  shock  wave  and 
confinement/water  interface  positions. 

Additional  assumptions  regarding  post  CJ  reacticn  are  then  necessary  to 
match  measured  and  predicted  positions.  Johnso’  and  co-workers  {3|  have 
applied  this  approach  to  ANFO,  alumiruzed  ANTO  and  several  other  ' 
commercial  products.  For  11%  aluminized  ANFO  the  results  could  be 
described  by  assuming  that  60%  of  the  AN  and  all  of  the  aluminium  react  by 
the  CJ  plane,  with  the  remaining  AN  reacting  within  a  few  microseconds. 
Similarly,  their  calculations  for  ANFO  indicated  that  55%  of  the  AN  remained 
inert  at  the  CJ  plane,  and  then  complete  reaction  occurred  within  a  few 
microseconds. 

However,  these  studies  employed  an  artificial  bum  model  to  advance  the 
detonation  front  in  the  hydrocode  simulations.  This  imposed  an  artificially 
large  radius  of  curvature  on  the  detonation,  front,  distorting  the  hydrodynamic 
flow  in  the  vicinity  of  the  explosive/confinement  interface.  Furthermore,  this 
procedure  imposed  an  artificially  short  reaction  zone  length,  and  hence  ignored 
the  energy  delivered  to  the  surroundings  by  expMsion  within  the  reaction  zone. 
By  performing  a  proper  two-dimensional  resolved  reactive  flow  simulation 
using  the  hydrocode  DYNA2D  modified  to  include  the  CPEX  equation  of  state 
and  reaction  kinetics,  Kennedy  (231  was  able  to  reproduce  the  ANFO  aquarium 
data  of  Johnson  et  al.  (31  with  no  a  priori'  assumptions.  The  predicted  CJ  zone 
length  was  16  mm,  with  the  o.dent  of  reaction  at  the  CJ  surface  being  95% 

(cf.  55%  from  Johnson  ef  al.)  on  the  charge  axis,  dropping  to  83%  at  the 
ANFO/confinement  boundary. 

Wannig>.r  et  al.  (24]  have  investigated  the  detonation  properties  of  a  variety  of 
ammonium  perchlorate  explosives  using  an  approach' sitiular  to  Johnson  el  ail. 
They  have  varied  the  fraction  of  AP  reacting  in  the  detonation  front  in  a  BKW 
calculation  until  agreement  is  obtained  with  experimental  detonation  velocity 
measurements.  They  have  then  combined  these  reaction  schemes  with  an 
underwater  hydrocode  to  calculate  relative  shock  and  bubble  energies. 

The  manner  in  which  we  have  measur^  the  reactioits  of  the  various 
components  of  PBXW-llS  to  die  reaction  rate  model  in  CPE.X  is  similv  to  the 
schemes  used  by  Murphy  et  al.  (251  and  Tarver  and  Green  (261  to  modd 
composite  explosives  and  propellants.  Murphy  et  al.  have  studied  a  family  of 
explosives  containing  HMX/AP/ZrH2/binder  using  a  three  term  ignition  and 
growth  model  The  first  term  in  the  expression  for  the  rate  of  reaction  is  the 
hot  spot  ignition  term  and  their  model  assumes  that  the  initial  hot  spot  fraction 
is  equal  to  the  initial  void  fraction.  The  second  term  describes  the  growth  of 
the  ignited  hot  spots;  previously  published  rates,  for  HMX  were  used  for  this. 


The  third  term  describes  the  rate  of  reaction  of  the  remaining  components. 
Tarver  and  Green  have  studied  a  propellant  composition  containing 
HMX/AP/Al/bindor  using  a  similar  scheme;  the  first  term  describes  hot  spot 
ignition,  the  second  the  giOwLh  rate  of  these  hot  spots  •;  ing  previously 
published  HMX  rates,  and  the  third  term  describ"?  tin-  ;;!■  bal  reaction  for  the 
rest  of  the  propellant  (AP/Al/binder).  The  unki. '  .-n  constants  in  the  reaction 
rate  terms  were  estimated  by  fitting  to  expterimental  data  from  embedded 
multiple  manganin  pressure  gauges  in  sustained  pulse  experiments  using  a  four 
inch  gas  gun.  Tarver  and  Green  then  used  this  calibrated  reaction  rate  model 
with  the  DYNA2D  hydrocode  to  estimate  the  failure  diameter  for  this  particular 
composition. 

Leiper  e-  il.  [81  have  also  used  the  CPEX  reaction  rate  moticl  in  the  DYNA2D 
code  to  model  an  air  sensitized  water-in-oil  emulsion  explosive  of  composition 
78.7%  AN,  16%  water  and  5.3%  oils  and  surfactants.  They  calculated  density 
profiles  and  other  flow  field  variables  and  found  good  agreement  with 
experiment.  '  ■ 

Interest  at  MRL  is  centred  on  the  cstiniation  of  relative  shock  and  bubble 
energy  for  PBXs  of  the  PBXW-115  type.  One  possible  approach  to  this  problem 
is  to  use  the  reaction  rate  parameters  obtained  from  CPEX,  calibrating  these  to 
experimental  data  on  the  unconfined  detonation  velocity-diameter  effect,  and 
then  using  this  reaction  scheme  in  a  one-dimensional  hydrocodc  to  simulate 
underwater  blast.  A  one-dimensional  hydrocode  for  detonation  calculations 
has  recently  been  developed  at  MRL  127]  and  is  currently  being  mod'fied  for 
application  to  problems  of  this  type  I28|.  Such  an  approach  however  requires 
the  availability  of  a  code  such  as  CPEX,  and  experimonta:  data  for  calibration 
•  purp  'ses.  It  should  be  poted  that  CPEX  took  10  nrian  years  to  develop  and, 
while  a  similar  code  could  probably  now  be  developed  in  less  time,  it  is 
unlikely  that  MRL  could  devote  the  necessary  resources  to  this. 

Another  approach  would  be  to  follow  the  method  of  analysis  of  Johnson  et  ai 
[3].  This  has  the  advantage  of  the  ready  a'-ailability  of  all  necessary  codes,  but 
would  require  the  establishment  of  facilities  to  en.iVle  aquarium  test  data  to  be 
obtained  for  PBXW-115.  It  should  be  noled  also  ij'at  Bdzil  and  co-workers 
(29,  30)  have  c'eveloped  a  slightly  different  analysis  of  diveigenf  flow  theory 
based  on  the  original  Wood-Ktrkwood  paper  aiiu  have  developed  a  Detonation 
Shock  Dyrtamics  (DSD)  code  to  model  flows  for  which  the  local  radius  of  . 
curvature  is  very  much  greater  than  the  reaction  zone  width  [31|. 

Whichever  approach  is  adopted,  it  is  evident  that  far  more  detailed 
experimental  characterization  of  the  particular  explosive  under  study  is 
required.  Hydrocodes  are  unable  to  make,  reliable  predictions  of  relative  shock 
and  bubble  energies  unless  considerable  experimental  data  is  available  to  allow 
accurate  parameterization  of  the  appropriate  reaction  rate  equations.  Once  this 
initial  information  has  been  obtained  modelling  techniques  can  then  be  used  to 
investigate  the  effect  of  small  changes  in  the  formulation  on  the  underwater 
performance  of  the  explosive.  These  techniques  will  decrease  the  fr^uirement 
for  expensive  field  testing  and  provide  efficient  and  relatively  fast  methods  for  ■ 
estimating  trends  in  underwater  performance. 
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